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ABSTRACT 

We analyse near-infrared fiST/NICMOS FllOW(J) and F160W(H) band pho- 
tometry of a sample of 27 i'-drop candidate z ~ 6 galaxies in the central region of 
the HST/ACS Ultra Deep Field (HUDF). The infrared colours of the 20 objects not 
affected by near neighbours are consistent with a high redshift interpretation. This 
suggests that the low redshift contamination of this i'-drop sample is smaller than 
that observed at brighter magnitudes where values of 10-40% have been reported. 
The J — H colours are consistent with a slope flat in /„ (fx oc A~ 2 ), as would be 
expected for an unreddened starburst. There is, however, evidence for a marginally 
bluer spectral slope (f\ oc A~ 2,2 ) which is perhaps indicative of an extremely young 
starburst (~ 10 Myr old) or a top heavy initial mass function and little dust. The 
low levels of contamination, median photometric redshift of z ~ 6.0 and blue spectral 
slope, inferred using the near-infrared data, supports the validity of the assumptions 
in our earlier work in estimating the star formation rates and, that the majority of 
the i-drop candidates galaxies lie at z ~ 6. 
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1 INTRODUCTION 

In recent years the observational horizon has expanded 
rapidly and radically for those observing distant galaxies. 
Large format red-sensitive detectors on wide field imaging 
instruments, the new generation of 8m class telescopes and 
the refurbished Hubble Space Telescope (HST), have pushed 
the limits to which we can routinely detect star-forming 
distant galaxies progressively from redshifts of one to 
beyond z = 6. At the highest redshifts currently accessible, 
narrow band emission lines searches using the Lyman-q line 
have moved on from redshifts of 4 llHu fc McMahorJll996t 
Fuiita et all 12003ft . to 5.7 (iHu. McMahon fc Cowie 
1999) and now reach to z ~ 6.5 IIKodaira et al. 



19991 iBremer et al. I 12004ft . spectroscopically confirmed, 



z ~ 6 


iStanwav. Bunker fc McMahon ll2003UStanwav et al. 


2004a 


Ibl iBouwens et al. 1 l2003allbl 12004: 


Dickinson et al. 


2004 


Yan. Windhorst & Cohen 112003ft . 



20031: IHu et al. I l2002j). Photometric redshifts 



Jjanzett^^^hi^fc^e^nand«z^o^J_^96j; JConnollv et al. I 
19971: I ? crnande z-Soto. Lanzetta fc Yahill 1 1999ft are now 
routinely used to process large datasets and identify high 
redshift candidates. An application of this method, the 
continuum based Ly man-break photometric technique 
fioneered at z ~ 3 by ( iGuhMhakiuita^ Tyson fc Maiewski I 



119901: ISteidel. Pettini fc Hamilton I Il995l) . has been ex- 
tended progressively to z~ 4 and z ~ 5 (e.g. ISteidel et al. I 



The 'classical' Lyman break technique used by Steidel 
et al. (1996) uses three filters, one redward of rest frame 
Lyman-a (A rcst > 1216A), a second in the spectral region be- 
tween rest-frame Lyman-a and the rest-frame Lyman limit 
(912A) and a third at A re st < 912A. At z~3 the technique 
relies on the ubiquitous step in the spectra of the stellar com- 
ponent of galaxies at 912A due to photospheric absorption, 
supplemented by optically-thick Lyman limit system absorp- 
tion caused by neutral hydrogen in the galaxy in question or 
in the intervening IGM. At higher redshifts the evolution in 
the Lyman-a forest absorption, particularly in the spectral 
region 912-1216A means that the effective break migrates 
redward to the Lyman-a region. Recently this search tech- 
nique has been extended to z ~ 6 by various authors and 
i'-drop samples have been used to constrain the star for- 
mation history of the universe (as derived from rest-frame 
UV luminous, starbursting stellar populations, e.g. Stanway, 



2 E. R. Stanway et al. 



Bunker & McMahon 2003; Giavalisco et al. 2004; Bouwens 
et al. 2004; Bunker et al. 2004). 

Any such estimation, however, is dependent on 
the characteristic properties of the i'- drop popula- 
tion. Although Steidel and coworker s dSteidel et al. I 
ll999l:ISteideL Pettini fc Adelberger1l200 ll) have successfully 
shown that a Lyman-break technique cleanly selects galax- 
ies at intermediate redshift (zsi3-4), the z'-drop varient 
of the technique lacks a second colour to constrain the red- 
shift distribution and is subject to contamination from both 
low redshift elliptical galaxies and cool dwarf stars. Until 
recently, these factors have been essentially unconstrained 
by observational data. 

The HST/ACS Ultra Deep Field (HUDF, Beckwith, 
Somerville and Stiavelli 2004), imaging an 11 arcmin 2 region 
of the sky to faint magnitudes in the F435W(B) , F606W (v) , 
F775W(i') and F850LP(z') bands, has now allowed a clean 
sample of i'-drop objects to be defined and constraints to b e 
placed upon its luminosity function dBunker et al. ll2004li . 
The Z/ST/NICMOS treasury programme, complimentary to 
the HUDF, has imaged the central region of the ACS field 
to faint magnitudes at wavelengths of 1.1 and 1.6 microns, 
allowing the infrared properties of the i'-drop population (in 
particular their rest-frame ultra-violet spectral slope, faint- 
end contamination and luminosity at 1500A) to be deter- 
mined. In this paper we discuss the infrared propert ies of 
the i'-drop population defined in lBunker et al. I i2004h . and 
their implications for the nature of these objects. 

We adopt a A-dominated, 'concordance' cosmology with 
tt A = 0.7, Q M = 0.3 and H = 70 /irokms -1 Mpc~\ All 
magnitudes in this paper are quoted in th e AB system 
llOke fc Gunn I Il983t and the iMadaul lll995li prescription, 
extended to z = 7, is used where necessary to estimate ab- 
sorption due to the intergalactic medium. Beyond z = 7 
absorption is assumed to decrease the transmitted flux by 
98% (the decrement observed in the flux of z > 6 quasars 
identified by the Sloan Digital Sky Survey, iFan et al. I2003T) 



2 HST OBSERVATIONS AND CATALOG 
CONSTRUCTION 

2.1 ifST/NICMOS Observations 

In this paper we utilise the publically released, pipeline- 
reduced images from the NICMOS UDF treasury pro- 
gramme 1 . These were taken with the NICMOS 3 camera 
on HST, in the FllOW('J') and FW0W( e H') filters, with 
exposures of 8 orbits per band per pointing giving an aver- 
age exposure time of 21500 seconds per pixel in each image. 
The transmission profiles of these filters are illustrated in 
figure Q and differ from those of standard ground-based J 
and H filters, tuned to atmospheric transmission windows. 
In particular the F110W filter is wide, extending down to 
8000A. 

NICMOS 3 uses a 256 x 256 pixel detector with a scale 
of 0'.'20 per pixel, giving a field of view of 51" x 51". The 
images discussed in this paper have been distortion corrected 

1 Programme GO-9803; PI Thompson; 144 orbit allocation; 
Available from http://www.stsci.edu/hst/udf 
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Figure 1. Normalised transmission profiles of the F775W(i') 
and F850LP(z') filters on the ACS instrument on HST and the 
FllOW(J) and F160W(H) filters of NICMOS. Note the signif- 
icant overlap between the J band and both lower wavelength 
filters. Overplotted on these are the profiles of a starburst galaxy 
placed at z ~ 6 (solid line, Kinney et al. 1996) and also a power- 
law profile wi th f\ oc A~ 20 (dashed line, see section 14.31 . The 
IMadaul Jl995fl prescription for intergalactic absorption at z = 6 
has been applied to both spectra. 

and drizzled onto a grid of 0'.'09 per pixel. The available 
NICMOS fields cover 2.3 x 2.3 arcmin giving a total area 
of 5.2 arcmin 2 (45% of that of the HST/ACS HUDF) to 
a depth of Jab = 27.73(3a) and Hab = 27.48(3ct) in a 2" 
diameter aperture. They comprise mosaics of 3 x 3 individual 
NICMOS pointings and the edge regions (with low exposure 
times) have been trimmed. Further deta ils on the use for 
NICM OS for deep surveys can be found in lThompson et al. I 
( 1999) which describes the results of 49 orbit J and H band 
observations of a sing le NICMOS pointing of the HDF-N 
jWillia.ms et al. B. 

The near-IR (J and H band) magnitudes of objects in 
these images are calculated using the photometric calibra- 
tion provided by the NICMOS UDF team, with photometric 
zeropoints given by Jab = 23.4034 and Hab = 23.2146 2 . 
The point spread function of unresolved objects (measured 
on the images) were f/.'4 FWHM (4 drizzled pixels) in the 
J-band and 0745 FWHM (5 drizzled pixels) in H. 

2.2 Catalog Construction 

Catalog construction from the imaging data was performed 
utili sing v2.3.2 of the SExtrac tor photometry software pack- 
age feertin fc Arnout3 [l996'l. For object identification, we 
demanded at least 5 adjacent pixels above a flux threshold 
of 2a per pixel. Catalogs were trained in the J band and 

2 Note that the photometric calibration of the NICMOS instru- 
ment appears to have changed significantly from that reported in 
earlier work (zeropoints of Jab= 22.7 and /f^s=22.8), probably 
due to the installation of a new cryostat. The older zeropoints 
have remained in use in the on-the-fly calibration software of the 
MAST archive. AB magnitudes may be converted to approximate 
Vega magnitudes usi ng the prescriptio n: Jab = •A'ega + 0.98, 
H AB = H Vega + 1.31 iKoornneef Il983tl 
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Sextractor was then used in two-image mode to evaluate 
the equivalent magnitude of each object in the H band. 

Objects in the resulting catalog were then matched to 
those in a second catalog of sources detected in the z' ACS 
HUDF image, with a match between two objects being as- 
signed if the astrometry agreed to better than (X'27 (3 pix- 
els on the drizzled J band image, one undrizzled NICMOS 
pixel) . The agreement between astrometric systems was gen- 
erally better than O'.'l although the discrepancy increases 
with distance from the field centre. 

Given the depth of the HUDF images, the narrower 
point spread function and smaller pixel scale ACS (PSF and 
plate scale are both C//05 for ACS compared to a plate scale 
of 0"2 and a PSF of 0'.'4 for NICMOS) can lead to difficulties 
with this process, since objects identified as single sources 
on the J and H band images may be resolveable into sev- 
eral distinct galaxies, or components of a single galaxy, on 
the optical imaging of the HUDF. As a result the z' flux 
corresponding to a single J-selected 'source' may be under- 
estimated. 

The SExtractor software performs source deblending 
based on a multi-threshold i ng alg orithm as described in 
detail in iBertin fc Arnout3 (Il996l) . This process involves 
rethresholding the flux of each set of connected pixels above 
the detection threshold into 30 levels spaced between that 
threshold and the peak flux. At each level, the source may 
be divided into two components if the fl ux in each compo- 
nent e xceeds a given threshold. As tests in lBertin fc Arnoutsl 
( 1996) indicate, such an approach may fail if there is no sad- 
dlepoint in the flux profile of the blended system (i.e. for 
objects seperated by less than two times the instrumental 
point spread function, or approximately 0.8" in this case). 
To verify this, we performed simulations in which a compact 
artificial galaxy of J = 25.0 was inserted at varying distances 
from real objects with 23 < J < 27 in the NICMOS UDF 
images. As expected, source magnitudes were successfully 
recovered for equal magnitude objects seperated by ~0'.'8, 
and for galaxies seperated by l'/2 even with a two magni- 
tude difference between real and artificial galaxies. On the 
other hand, the software could not deblend objects closer 
than these seperations. Hence the photometry of i'-drops 
with brighter neighbours within 1" (indicated with an as- 
terisk in table is not used in the analysis in sections 2] & 

m 

We use deblended "total" magnitudes, rather than fixed 
aperture magnitudes in the J and H bands. SExtractor cal- 
culates the "total magnitude" by computing an elliptical 
aperture whose ellipticity and elongation are defined by the 
second order moments of the object's isophotal light distri- 
bution. The flux is measured in an elliptical aperture with 
2.5 times the elongation calculated above. Hence the mea- 
sured flux is independent of an assumed radial light profile. 

The use of different photometric apertures in different 
observation bands may lead to low surface brightness flux 
being missed in one or more bands. We attempted to quan- 
tify this by comparing the z-J colours obtained using "total" 
magnitudes with those obtained using fixed circular aper- 
tures, and corrected for missing flux assuming a compact 
radial light profile as appropriate in the z' band. We found 
that the resulting colours agree to within approximately 0.2 
magnitudes. Given the consistent colours obtained, we de- 
cided to adopt the "total" magnitudes rather than assuming 



a radial light profile for this source population. However, we 
note that without knowledge of the accurate radial light pro- 
file some uncertainty remains in the fraction of flux emitted 
below the detection threshold in each band. 



3 NEAR-IR PROPERTIES OF Z ~ 6 
CANDIDATE OBJECTS 

3.1 Near-IR properties of bright I-drops in the 
HUDF 

In a previous paper dBunker et al. Il2004f) . we reported a 
catalogue of 54 z ~ 6 galaxy candidates, including one con- 
firmed z = 5.83 galaxy (Bunker #20104, SBM03#1, SiD002, 
GLARE 1042 in Bunker et al. 2004, Stanway et al. 2004a, 
Dickinson et al. 2004 and Stanway et al. 2004b). Each i'- 
drop object in this selection satisfied (i' — z')ab > 1.3 and 
z' AB < 28.5 in the HST/ACS images of the HUDF field. Of 
this catalog, 27 objects lie in region surveyed by the NIC- 
MOS UDF (consistent with the 45% of the ACS HUDF area 
covered by this near-IR imaging.) 

The near-IR photometry of these high redshift candi- 
dates is presented in table Q with 3<r limits on the magni- 
tudes indicated where appropriate. The photometry of seven 
of the objects is likely to be affected by bright near neigh- 
bour objects at lower redshift. 

Interestingly several multicomponent systems of galax- 
ies are observed, perhaps implying that we are observing 
merger-triggered star formation. As illustrated by figure |K| 
three i'-drop galaxies, 24228, 23972 and 24123, form a close 
system (spanning m 1") and hence have not been resolved 
as individual objects in the near-IR (where each pair sepa- 
rated by just 3 undrizzled pixels). As a result, the average 
colour of this galaxy group, Group 1, has been calculated 
by summing the optical fluxes of the three components. We 
note that the majority of near-IR emission is to the north 
of the group centre, suggesting that object 23972 may be 
fainter than its neighbours in the near-IR bands. 

Similarly, the galaxy pairs 44154 & 44194 (Group 2) 
and 42929 & 42806 (Group 3) are each treated as a sin- 
gle system, their optical flux summed and compared to the 
blended near-IR flux of each system. Hence the reported flux 
for each of these groups of z'-drop galaxies is a system total 
rather than that of any particular galaxy. 

3.2 Near-IR properties of fainter I-drops in the 
HUDF 

Fourteen candidate high redshift galaxies identified by 
Bunker et al. (2004) lie within the NICMOS field yet are 
undetected to 3a in the J and if-band images. Of these, 
nine are isolated, lying more than an arcsecond away from 
their nearest neighbours. In order to investigate the proper- 
ties of this fainter sub-catalogue of i'-drop objects, sections 
of the NICMOS images, centered on the z'-band coordinates 
of the 9 isolated candidate objects, were stacked. The resul- 
tant images (shown in figure^! were then analysed with the 
SExtractor software as described in section |2~2*1 

The background noise in the stacked images is a factor 
of y9 lower than in the original NICMOS images, allowing 
the average properties of these candidates to be probed to 
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Dec (J2000) 








z 


AB 




Jab 


Hab 


(z' - J) AB 


{J - 


H)ab 






Group l 1 


03 


32 


34. 


.29 


-27 


47 


52. 


.8 


28.06 ± 0.15 


26. 


.11 


± 





.05 


26.23 ± 0.10 


26.09 ± 0.11 


0.18 ± 0.11 


0.14 


± 


0.15 


5. 


qo+0.24 
ao -0.26 


Group 2 2 


03 


32 


36. 


.46 


-27 


46 


41. 


5 


28.72 ± 0.22 


26. 


.36 


± 


0. 


.05 


26.66 ± 0.11 


25.82 ± 0.07 


-0.33 ± 0.12 


0.84 


± 


0.13 


5. 


' J -0.12 


Group 3 3 


03 


32 


37. 


.46 


-27 


46 


32. 


.8 


29.74 ± 0.40 


27. 


.27 


± 





.09 


26.53 ± 0.13 


26.29 ± 0.13 


0.94 ± 0.16 


0.24 


± 


0.18 


6. 


fi c:+0.35 
D0 -0.11 


20104 


03 


32 


40. 


.01 


-27 


48 


15. 





26.99 ± 0.04 


25. 


.35 


± 





.02 


25.54 ± 0.04 


25.51 ± 0.05 


-0.19 ± 0.04 


0.03 


± 


0.06 


5. 


fio +0.01 
SZ -0.01 


23516* 


03 


32 


34. 


.55 


-27 


47 


56. 





28.57 ± 0.10 


27. 


.04 


± 





.05 


26.82 ± 0.13 


26.85 ± 0.16 


0.22 ± 0.13 


-0.03 


± 


0.21 


5. 


70+O.O8 
'°-0.06 


25941 


03 


32 


33. 


.43 


-27 


47 


14. 


9 


29.30 ± 0.18 


27. 


.32 


± 





.06 


27.52 ± 0.16 


27.24 ± 0.15 


-0.20 ± 0.17 


0.28 


± 


0.22 


5. 


Ql+0.09 
a± -0.07 


26091* 


03 


32 


41. 


.57 


-27 


47 


14. 


2 


29.74 ± 0.25 


27. 


.38 


± 





.06 


26.44 ± 0.12 


26.33 ± 0.13 


0.94 ± 0.13 


0.11 


± 


0.18 


6. 


fi7 +0.15 
"'-0.06 


24458+ 


03 


32 


38. 


.28 


-27 


47 


51. 


3 


29.11 ± 0.15 


27. 


.51 


± 





.07 


>27.7 


>27.4 


<0.19 








5. 


oi +0.09 
ol -0.06 
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03 


32 


38. 


.96 


-27 


46 


00. 


5 


29.77 ± 0.26 


27. 


.74 


± 





.08 


26.36 ± 0.10 


25.35 ± 0.06 


1.34 ± 0.13 


1.01 


± 


0.12 


6. 


on+0.05 
-0.05 
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03 


32 


35. 


.06 


-27 


47 


40. 


2 


>30.4 


27. 


.83 


± 





.08 


27.32 ± 0.16 


27.30 ± 0.20 


0.51 ± 0.18 


0.02 


± 


0.26 


6. 


o,+0.30 
°*-0.15 


14751+ 


03 


32 


40. 


.92 


-27 


48 


14. 


.8 


29.39 ± 0.17 


27. 


.87 


± 





.09 


>27.7 


>27.4 


<0.14 








5. 


7 o+0.13 
'°-0.07 


35084* 


03 


32 


41 


.70 


-27 


47 


11. 


6 


29.86 ± 0.28 


27. 


.92 


± 





.09 


>27.7 


>27.4 


<0.22 








5. 


qo+0.15 
ao -0.17 


46503+ 


03 


32 


38. 


.55 


-27 


46 


17. 


5 


29.43 ± 0.20 


27. 


.94 


± 





.09 


>27.7 


>27.4 


<0.24 
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lyyoo 
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"1 A 


a 
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97 
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no 
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0. 


7O+0.15 
■ '°-0.09 


2111lt 


03 


32 


42. 


.60 


-27 


48 


08. 


9 


29.69 ± 0.24 


28. 


.02 


± 





.10 


>27.7 


>27.4 


<0.32 








5. 


on+0.23 
—0 11 


46223+ 


03 


32 


39. 


.87 


-27 


46 


19. 


1 


>30.4 


28. 


.03 


± 





.10 


>27.7 


>27.4 


<0.33 








6. 


U °-0.32 


22138+ 


03 


32 


42. 


.80 


-27 


48 


03. 


3 


>30.4 


28. 


.03 


± 





.10 


>27.7 


>27.4 


<0.33 








6. 


n o+0.23 
U °-0.32 


46234 


03 


32 


39. 


.86 


-27 


16 


19. 


1 


>30.4 


28. 


.05 


± 





.10 


>27.7 


>27.4 


<0.35 








6. 


08 +0.24 
u -0.33 


12988* 


03 


32 


38. 


.49 


-27 


48 


57. 


.8 


>30.4 


28 


.11 


± 





.11 


>27.7 


>27.4 


<0.41 








6. 


ns +0.25 
US -0.35 


24733+ 


03 


32 


36 


.62 


-27 


17 


50. 





>30.4 


28 


.15 


± 





.11 


>27.7 


>27.4 


<0.45 








6. 


ns +0.25 
U °-0.36 


21530+ 


03 


32 


35 


.08 


-27 


48 


06. 


.8 


30.24 ± 0.39 


28 


.21 


± 





.12 


>27.7 


>27.4 


<0.51 








5. 


og+0.21 
a -0.32 


35271+ 


03 


32 


38 


.79 


-27 


47 


10. 


9 


29.77 ± 0.26 


28 


.41 


± 





.14 


>27.7 


>27.4 


<0.74 








5. 


7n +0.25 
,u -0.13 


22832* 


03 


32 


39 


.40 


-27 


47 


59. 


.4 


>30.4 


28 


.50 


± 





.15 


>27.7 


>27.4 


<0.80 








6. 


n o+0.34 
U8 -0.50 


Stack (mean) + 
















30.14 ± 0.17 


28. 


.13 


± 





.05 


28.56 ± 0.16 


28.45 ± 0.16 


-0.43 ± 0.17 


0.11 


± 


0.23 







* Photometry may be affected by near neighbour objects. 

+ Object has contributed to the stack of undetected objects, the mean properties of which are given in the last line. 

1 2 3 J and H band fluxes of these groups of objects are blended (see text for details). Photometry is given for the group as a whole. 

Group 1 comprises objects 24228, 24123 & 23972. Group 2 comprises objects 42929 & 42806. Group 3 comprises objects 44154 & 44194. 



Table 1. Photometry of lBunker et al.1 l2004fl candidate objects lying within the NICMOS UDF Field. Near-IR magnit udes given are 
total m agnitudes, i' and z' magnitudes were measured in 0'.'5 apertures and corrected to total magnitudes as described in lBunker et al.1 
(2004), except for those cases where the i'-drop galaxies form small groups. In these cases the total magnitude of the group is reported. 
U pper limits (3a , measu red in a 2" diameter aperture) are given for the J and H bands where appropriate. ID numbers are those given 
in lBunker et alH l2004Tl . Photometric redshifts were calculated as described in section 14.21 and the 99% (3cr) confidence intervals are 
shown. 



a depth over a magnitude fainter in each band. Their mean 
photometry is shown in table Clearly, this analysis as- 
sumes that the 9 i'-drop z ~ 6 candidates considered are 
similar in colour and that the mean colour is not heavily 
skewed by a single source. However, the objects analysed 
are similar in z flux (spanning less than 1 magnitude) and 
have near-IR limits in the NICMOS images which rule out 
identification as low redshift elliptical galaxies. The mean 
colours of these faint galaxies are bluer than those for which 
secure NICMOS detections are available, suggesting a steep 
mean rest-frame UV spectral slope (see section I47JI . 



4 DISCUSSION 

4.1 Faint Contaminants in i'-Drop Selection 

The very red, old stellar populations of elliptical galaxies at 
around 2 = 1 — 2 (Extremely Red Objects, or EROs) allow 
them to contaminate any sample of objects selected using a 
simple i — z colour criterion (as illustrated by figure [5] sec 
also figure 2 of Stanway, Bunker & McMahon 2003). The 



degree of sample contamination has long been uncertain, af- 
fecting the analysis of population statistics at high redshift. 

Utilising the v0.5 release of data from the deep 
Great Observatories Origins Deep Survey (GOODS), we 
have previously defined a sample of 20 i'-drop ob- 
jects with z' < 25.6 in the 300 arcmin 2 GOODS 
fields iStanwav. Bunker fc McMahon ll2003t IStanwav et al. I 
l2004al) . Of these 2 (10%) are probable EROs and a fur- 
ther six (30%) cool Gala ctic stars giving a cont amination 
at the bright end of 40%. iDickinson et al. I 1120041) . working 
with somewhat deeper GOODS data and to a lower signal- 
to-noise limit reported a contamination fraction of their i'- 
drop sample that increased at fainter magnitudes, reaching 
as high as ~45% for resolved galaxies (c.f. 10% for our bright 
sample), although they report that their catalogue was af- 
fected by a number of spu rious sources no t expe cted to be 
a problem in the HUDF. iBouwens et al. I <l2004l) . also us- 
ing GOODS (vl.0) and with a similar selection, reported a 
contamination of their i'-drop sample by lower redshift el- 
lipticals of 11% based on the presence of flux in the B or V 
bands (unlikely for galaxies at z ~ 6). 

The availability of NICMOS imaging in the central por- 
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Figure 2. z' — J v i' — z' Co l our-co lour diagrams. The evolutionary tracks followed by a starburst spectrum iKinnev et al, Il996t) and the 
IColeman. Wu fc Weedman I ll98dl) empirical galaxy templates are plotted as lines (as in legend of of the right hand figure) . The colours 
of the z ~ 6 candidate galaxies of Bunker et al. (2004) are indicated with solid symbols. Colours of the galaxy groups as described in the 
text are indicated by diamonds. Objects for which only 3cr limits are available are shown separately on the right hand figure for clarity. 
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Figure 3. J — H v z' — J Colour-colou r diagr ams. The evolution- 
ary tra cks followed b y a[Kirmev^^al A ^19961) starburst spectrum 
and the lColeman. Wu^^Ve^dmanTTl^80^ eiiiDirical galaxy tem- 
plates are plotted as lines (as in figure In addition tracks for 
a synthetic spectra with power law slopes of /3 = —2.0 and —2.5 
(fx = \&) are sho wn (long dashes). The colours of the z ~ 6 can- 
didate galaxies of iBunker et al.1 120041) are indicated with solid 
symbols. The colours of i'-drop groups are indicated with dia- 
mond symbols. 
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Figure 4. Stacked images of the nine objects undetected in the 

NICMOS J and H band images which are more than 1" from their 

nearest neighbour. The combined i'-drop sample is well detected 

in both near-4R wavebands. The regions shown are 3" by 3" in 

size 
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tion of the HUDF field allows lower redshift galaxies to be 
cleanly excluded from our sample, as illustrated in figures|31 
andU Of the 27 i'-drop objects for which NICMOS infrared 
data is available, only one (49117) has z — J colours which 
may be consistent with those of z ~ 2 elliptical galaxies. 

Object 49117 is likely tp be subject to contamination 
from two neighbouring objects: one star, bright in the J 
band (Jab = 22) and separated from the candidate by 2'.'3, 
and a second galaxy of comparable J magnitude to the can- 
didate and separated from it by only 0'.'6 (< 2 x FWHM, 
see figure [5J, unresolved from it. None of the objects unde- 
tected in J and H have colour limits consistent with lying at 
low redshift. This allows us to place an upper limit on the 
low redshift galaxy contamination in this i'-drop selected 
catalog of at most 1 objects in 27 (< 4%), or more likely no 
z ~ 2 elliptical galaxies in this sample. 

All these objects are also resolved in the z'-band, the 
only unresolved i'-d rop object reported in this field by 
IBunker et al. I <l2004l) having been excluded for having sig- 
nificant V-band flux. When combined with the near-IR evi- 
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Figure 5. Postage stamp images from the ACS and NICMOS for the first 10 objects listed in TableQ] The regions shown are 3" by 3" in 
size and orientated with north upwards. As can be seen, object 49117 is contaminated in several bands by flux from its near neighbours 
and a diffraction spike, while the groups of objects 24123/24228/23972, 42929/42806 and 44154/44194 are each blended into one source 
in the near-IR bands (labelled 1, 2 and 3 respectively in table HI, 
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dence from the NICMOS UDF exposure, this suggests that 
the contamination of this i'-drop sample is less than 5% - 
significantly lower than has been found by previous i'-drop 
studies at brighter magnitudes, suggesting that we may be 
seeing through the Galactic disk at these faint magnitudes. 

4.2 Redshift Discrimination 

The i'-drop method is, in principle, sensitive to any galaxy 
bright in the rest-frame ultraviolet and lying at redshifts be- 
tween z — 5.6 and z = 7 (where the Lyman-break shifts out 
of the z filter). However, the sharp fall off in filter trans- 
mission beyond 8500A biases the selection towards the low 
redshift end and renders unlikely the detection of galaxies 
beyond z w 6.5. 

The availability of extremely deep imaging across a wide 
wavelength range, allows the photometric redshift method to 
be used with go od results. The publically av ailable software 
package hyperz feolzonella. Miralles fc Pellct 2000) was used 
to calculate the photometric redshifts of each object in ta- 
ble Q given their magnitude, or magnitude limit, in the v, 
i 1 , z', J and H bands (spanning A rest = 4000 - 18000 A). 
This package performs x 2 fitting of the object flux to a set 
of predefined template spectra, selecting both the best fit- 
ting template and the best fitting redshift simultaneously. 
IColemam Wu fc Weedman I Jl98Ctl E, Im, and Sbc tem- 
plates were used, together with power law spectra of the 
form fx oc A" with =-1.0, -1.5, -2.0, -2.2, -2.5 (where 
— —2.2 is suggested by the infrared colours of the i'-drop 
sample, see section FOl . 

Trial redshifts were varied in the range < z p hot < 7 at 
intervals of Az — 0.05. All sources were found to lie above 
£phot > 5.6 using this method. To improve redshift discrim- 
ination the range was further narrowed to 5 < z p hot < 7 for 
the remaining objects, with a redshift interval Az = 0.02, re- 
sulting in the photometric redshifts presented in tableland 
illustrated in figure |S| All objects were best fit by power law 
spectra with slopes > —2.0 with the exception of objects 
49117 (CWW E). 

One i'-drop in this field has a confir med spectroscopic 
redshift, object 20104 at z svec = 5.83 iDickinson et al. I 
12004 IStanwav et al. Il2004albh . Its photometric redshift is 
consistent with this (z p hot = 5.82 ±0.01). The mean z p hot of 
the total i'-drop sample is 6.02 ±0.30. The mean z p hot of the 
i'-drop sample using only the objects(16) that have photom- 
etry unaffected by blending is 5.96 ± 0.22 and 6.16 ± 0.43 
for those that may be effected by blending with compan- 
ion objects(7). Hen ce a central assumed redshift of 6.0 
( Bun ker et al. 1 12004) is in good agreement with the data. 

4.3 Effect of Spectral Slope and Line 
Contamination 

Interestingly, those objects for which measured colours are 
available indicate an average spectral slope steeper than = 
—2.0, as do the mean colours of fainter objects detected in 
the stacked image. The evidence for this trend is clear in 
figure[3|as the J — H colour has a longer baseline than z' — J 
and hence is more sensitive to spectral slope. Although the 
error on each colour is large, the colours of most of those 
objects for which photometry is available lie between the 
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Figure 6. Photometric redshifts for the i'-drops in table HI The 
dotted lines are for the sample without neighbours and the dashed 
line is for objects with photometry that may be effected by 
neighbours. These redshifts were de rived using the hyperz tem- 
plate fi tting package and a mixture o flColeman. Wu fc Weedman I 
( 1980) templates and power law spectra with slopes in the range 
1.0 </3 < 2.5. 



redshift tracks predicted for objects with spectral slopes of 
= -2.5 and = -2.0. 

Strong line emission in the J-band might explain these 
observed colours. Narrow-band selected Lyman-a emitters 
at z ~ 5.7 are known to have strong emission lines, 
with typical rest-frame equivalent widths (EWs) of 50-150A 
iAiiki et al. 1120031 : iHu et al. Il2004|) and in some cases much 
higher (H^est > 200A e.g. Malhotra & Rhoads 2002, Rhoads 
et al. 2003). In the observed frame at z ~ 6, these correspond 
to equivalent widths of ~1000A. The F1WW( J)-band filter 
on NICMOS is broad with a FWHM of 5880A. Hence a line 
emitter with an observed equivalent width of 1000A could 
contribute up to 17% of the flux in this band or 0.2 mag. 
This, in fact, underestimates the potential contribution of 
line flux to the J-band magnitude since line blanketing be- 
low the Lyman break will suppress the J-band flux for ob- 
jects at redshifts greater than 5.6. 

The observed colours are thus consistent with a red- 
der intrinsic rest-frame UV spectral slope (due to stellar 
continuum emission) of — —2.0, if every object has a 
powerful emission line. The single i'-drop object in this 
field with a spectroscopically confirmed emission line, how- 
ever, has a much lower rest-EW (VKrost ~ 30A, contribut- 
ing <5% of the z flux, Stanway et al. 2004a, Dickinson et 
al. 2004), as do several other i'-drop objects which have 
now been confirmed s pectroscopically l|Bunker et al. ll2003l : 
IStanwav et al. Il2004bh . 

An alternative explanation for the observed near-IR 
colours of this sample may be evolution in the intrinsic spec- 
tral slope of starbursts from that observed at lower redshifts. 
The observed rest-frame UV spectral slope in star form- 
ing galaxies is influenced by two factors: the intrinsic spec- 
tral energy distribution of the stellar population (shaped 
by the initial mass function (IMF), and star formation his- 
tory) and the reddening due to dust extinction. The Lyman 
break population at z — 3 has a typical observed spectral 
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slope of = -1.5 ± 0.4 llAdelberger fc Steidell Eoocl) . as 
expected for dust-reddened starbu rst galaxies with an in - 
trinsic spectral slope of j3 = —2.0. ILeitherer et, al. I |l999), 
using their Starburst99 stellar synthesis code, find that star- 
bursts less than ~ 200 Myr in age are well approximated by 
the power law (f\ oc A' 3 ) model. They predict an intrinsic 
spectral slope at 1500A of f3 = -2.5 for young (1-10 Myr), 
low metallicity starbursts (Leitherer et al. 1999, figures 71 
& 72), independent of evolution and initial mass function 
(IMF) eff e cts (w here IMFs considered include the classical 
ISalpeter I Jl955f) IMF, a Salpeter IMF truncated at 30M Q 
and a steeper IMF with an enhanced proportion of low mass 
stars). In these young systems, hot OB type stars dominate 
the rest frame ultraviolet continuum. The predicted j3 falls 
to -2.0 on an IMF dependent timescale of the order of 100 
Myr and falls more slowly with increasing metallicity. 

These models do not account for the contribution to 
the UV spectral slope due to dust reddening. The typical 
colour excess fro m the stellar continu um of Lyman break 
galaxies at z ~ 3 JSteidel et al. Ill999b is E S (B - V)=0.16. 
ICalzetti et al. I i2000T) . studying a large sample of local 
starburst galaxies, found a relation between colour excess, 
E(B — V) and observed spectral slope, /3, of E(B — V) = 
0.52(/3 — 0o) where /3o is the intrinsic spectral slope, although 
the scatter on this relation is considerable (approx 0.5 mag). 
If the dust properties of the very young starbursts at z ~ 6 
are similar to those of Lyman break galaxies at lower red- 
shift, this would imply an intrinsic spectral slope of -2.9 for 
objects with /3 = —2.2 (typical for our sample, see sections 
I4.2l and l4.31 . outside the range of slopes predicted from star- 
burst models. This suggests that, if the J — H colours - and 
photometric redshifts - are good indicators of spectral slope, 
this i'-drop sample is less affected by dust reddening than 
Lyman break galaxies at lower redshift. 

Given the fraction of i'-drops that are observed to be 
in close pairs which may be interacting systems, these near- 
IR observations suggest that we may be observing young 
starbursts triggered by an epoch of merger activity. 



5 IMPLICATIONS FOR STAR FORMATION 
HISTORY 

The linear relation between star formation rate 
SFR) and UV luminosity density at A rcst = 1500 A 
Madau. Pozzetti fc Dickinson Ill998t) has been widely used 
to study the volume-averaged star formation history of 
the universe (e.g. Madau et al. 1996, Steidel et al. 1999). 
The broad-band photometry of a source may be used 
to determine its luminosity density, if the redshift and 
appropriate fc-correction (from the effective wavelength of 
the filter to 1500A) is known. Clearly the fc-correction for 
UV luminous galaxies depend on their spectral slope, /3. 

Galaxies selected by the i'-drop technique present sev- 
eral difficulties in this regard. Firstly, with no strong 
constraint from bands at longer wavelength than z 
(A e ff=9200A, A rcst = 1300A at z = 6), the spectral slope 
of these objects has been unknown and a value of j3 = — 2.0 
has been assumed (e.g. Stanway, Bunker & McMahon 2003) 
or a lower redshift Lyman-break galaxy has been adopted as 
a template (e.g. Giavalisco et al. 2004). Furthermore, since 
the redshift distribution of objects is much broader than that 



of galaxies selected using the more traditional two-colour 
Lyman-break approach, luminosity densities have been cal- 
culated by placing the candidate objects at some model- 
dependent redshift, representative of the sample as a whole. 
For i'-drop samples t his representative re dshift has var- 
ied between z = 6.0 feunker et al. 1 120041 and z = 5.74 
iGiavalisco et al. Il2004l) depending on details of models and 
the selection criteria. 

Clearly, if the majority of objects lie at redshifts much 
higher or much lower than has been estimated, the star for- 
mation density calculated in earlier works will be unreli- 
able. While the near-IR data now available for i'-drop can- 
didates in the NICMOS UDF do not place a tight constraint 
on the redshift of individual objects, their properties (dis- 
cussed in sections 14.21 and I4.3|l are reassuringly consistent 
with those p redicted in simulations for the redshift range 
z = 5.6-6.5 1st anw av et al. 1120 04b: Giavali sco et al. 1 12004; 
iBouwens et al. Il200^) ~with a median colour consistent with 
lying at z ~ 6 as expected. However, as discussed in sec- 
tion there is some evidence for a slightly steeper spec- 
tral slope than the assumed value of j3 = —2.0. A slope of 
(3 < —2.0 would imply that the fc-corrections hitherto ap- 
plied to the data have led to a slight overestimation of the 
star-formation density. 

Figure compares the mean SFRs derived from the 
z' and J band filters, centered at A rcst = 1400A and 
1600A at z = 6.0 respectively, for the subsample of ob- 
jects isolated from lower redshift sources and for which a 
secure J band detection was made. The l umin osity- SFR re- 
l ation of iMadaii^ Pozzetti fc Dickinson I l)l998f) (which uses 
a lSaloeter I l)l955f) initial mass function) is used (Luv — 8 x 
10 27 (SFR/M yr" 1 ) ergs s" 1 Hz -1 ). Object s were assumed 
to lie at the mean redshift of z = 6.0 (see iBunker et al. I 
l2004h . The total flux of these seven objects in each wave- 
band was summed before conversion to star formation rates 
using the spectral slope specified. If the value of /3 used to 
calculate the fc-correction accurately represents the popula- 
tion mean slope, then the star formation rates derived from 
any two bands should converge. 

The photometric errors at these faint magnitudes, 
and the consequent uncertainties in the SFRs, are large. 
Nonetheless, the slope which gives the best agreement be- 
tween z and J-band derived star formation rates in the 
i'-drop sample (/3 = —2.12 ± 0.15) does suggest a slightly 
steeper slope than /3 = —2.0, as expected given the J — H 
colour distribution. In practise, varying the spectral slope 
between f3 — —2.0 and /3 = —2.2 changes the calculated star 
formation rate by only ss 6%. The derived value of (3 is also 
consistent with a slope flat in /„ (f3 — —2.0) used in ear- 
lier computations of star formation rates (e.g. Bunker et al. 
2004). 

The infrared properties of this i'-drop sample do sug- 
gest further caveats to be applied when calculating the star 
formation density at z ~ 6. If, i'-band drops are indeed 
young starbursts, then the iMadau. Pozzetti fc Dickinson I 
(1998) luminosity-star formation rate relation may be 
underestimating the star formation density. Stellar pop- 
ulations <10 Myr old y ield SFRs more than 50% 
higher than given by the IMadau. Pozzetti fc Dickinson I 
(1998) cali bration when modelled using stellar synthesis 
techniques dLeit herer fc Heckman I Il995l : iKennicutt I Il99sl : 
ILeitherer et al. Ill999l) ~ 
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Figure 7. The ratio of SFRs derived from J and z' band fluxes 
for those objects with secure J band fluxes (objects 20104, 25941, 
27270 & groups 1 and 3) for different spectral slopes. Solid regions 
indicate the lcr uncertainty on star formation rates derived from 
z mag nitudes (A rost = 1400A), line-filled regions indicate star 
formation rates from the J-band (observed at A re st = 1600A). 
In each case the 1500A luminosity density has been corrected 
using the spectral slope indicated and assuming the objects lie at 
z = 6.0. 



On the other hand, a top-heavy IMF rich in massive 
stars has been hypothesised at very high redshifts (e.g. 
z > 10 - 15, Clarke & Bromm 2003) and such an IMF 
could explain the steep UV spectral slope observed. If so, 
this would also alter the luminosity- SFR calibration which 
extrapolates the formation of high mass stars to the low 
mass regime. A top heavy IMF would reduce the star for- 
mation rate inferred from a given 1500A luminosity. 



6 CONCLUSIONS 

(i) The near-infrared properties of a faint sample of i! - 
drop objects (defined by Bunker et al. 2004) has been inves- 
tigated, using the deep FllOW 7 and F1Q0W images of the 
HST/ACS Ultra Deep Field. 

(ii) The infrared colours of all objects not suffering from 
contamination by near neighbours are consistent with a high 
redshift interpretation. This suggests that the contamina- 
tion of this i'-drop sample is smaller than that observed at 
brighter magnitudes. 

(iii) The spectral slopes inferred from near infrared 
colours are consistent with f3 — —2.0 although there is 
evidence for a marginally steeper spectral slope of /3 = 
-2.2 ±0.2. 

(iv) These steep spectral slopes suggest that the dust ex- 
tinction of the sample is small and that the i'-drop popula- 
tion may comprise galaxies with young starbursts - possibly 
triggered by mergers (several close systems of galaxies are 
observed in this sample). 

(v) The low levels of contamination, and the steep spec- 
tral slope, inferred from near-infrared data supports the va- 
lidity of our previous star formation density estimates based 
on i'-drop samples. 
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